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Abstract 


The initial condition of high-mass star formation is a complex area of study because of the high densities 
(ny, > 10° cm ?) and low temperatures (Tgust < 18 K) involved. Under such conditions, many molecules become 
depleted from the gas phase by freezing out onto dust grains. However, the N-bearing and deuterated species could 
remain gaseous under these extreme conditions, suggesting that they may serve as ideal tracers. In this paper, using 
the Plateau de Bure Interferometer and Very Large Array observations at 1.3 mm, 3.5 mm, and 1.3cm, we 
investigate the possible habitats for NH3, NH2D, HCN, HCPN, SO, and CO in eight massive precluster and 
protocluster clumps G18.17, G18.21, G23.97N, G23.98, G23.44, G23.97S, G25.38, and G25.71. We found that 
the NH; cores are in good agreement with the 3.5 mm peak emission, but the NH; is much more extended than the 
3.5 mm emission structure. The SO distributions agree well with the 3.5 mm peaks for the evolved star formation 
stage, but we did not detect any SO emission in the four earliest star formation sources. C'?O is a poor tracer in 
conditions of the cold (<18 K) and dense (Z;10* cm ?) cores, e.g., the prestellar cores. We also found that the 
NH2D cores are mainly located in the temperature range of 13.0-20.0 K, and the NH5D lines may be strongly 
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1. Introduction 


High-mass star formation with larger than 8 Mẹ dominates 
the environment and evolution of the Galaxy. However, there is 
still a lack of knowledge about how high-mass stars form in 
their early stages (e.g., Pillai et al. 2007, 2011, 2012; Motte 
et al. 2018; Zhang et al. 2019, 2020). Particularly, the transition 
from a prestellar core to a protostar has not been well studied 
observationally due to the short lifetime of the prestellar stage, 
the difficulty in determining the relative evolutionary stage of 
the objects, the small spatial scales involved, and the lack of 
strong molecular tracers of cores at high densities (Friesen 
et al. 2014). At the high densities (ny, > 106 cm? ) and low 
temperatures (Tau. < 18 K) characteristic of the interiors of 
highly evolved star-forming cores, many molecules deplete 
from the gas phase by freezing onto dust grains (Walmsley 
et al. 2004; Bergin & Tafalla 2007), even the dense tracer co 
(Zhang et al. 2017). However, some N-bearing species are 
good tracers of dense gas and cold conditions toward initial star 
formation. In addition, the deuterium fraction of the remaining 
gas phase species increases dramatically above the cosmic 
abundance ratio ([D/H] ~ 1.5 x 1075; Oliveira et al. 2003) due 
to the increased production of HD * by the reaction of H} with 
HD in places where CO is depleted (Roberts & Millar 2000). 


The elements C, N, O, and S, produced by nucleosynthesis 
in stars, are the important components in the early star 
formation of the Galaxy (e.g., Wilson & Rood 1994; Wielen & 
Wilson 1997; Ikeda et al. 2002). The abundances of C, N, and 
O may contain important information not only about the history 
of nucleosynthesis, but also about the earliest stages of star 
formation. Understanding their precise content in local 
interstellar clouds is of fundamental importance in astrophysics 
and astrochemistry (e.g., Gerner 2014; Gerner et al. 2015). The 
elements have been widely observed in many different 
observations, such as in the form of CO, SO, HCN, and NH3 
(e.g., Zhang et al. 2014, 2017). However, their isotopic species, 
such as SO, HPÜCN, HCPN, and NH;D (e.g. Daniel 
et al. 2016; Booth et al. 2019; Li et al. 2024) in the early 
stage of massive star formation, are rarely studied in detail. 

Investigating different molecular tracers in different star 
formation environments is helpful to tell us how to select and 
use an appropriate tracer to study a particular star formation 
condition. Thus, in this paper we mainly investigate the 
possible habitats for NH3, NH5D, H?CN, HCPN, SO, and 
C'*O associated with the early stage of high-mass star 
formation. In addition, in Section 2 we present the associated 
observations and data reduction. In Section 3 we present the 
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observational results and analysis for the 3.5mm continuum 
emission, NH3, NH;D, H^CN, HCN, SO, and C!80 spectra. 
In Section 4 we discuss possible habitats for NH5D and the 
other N-bearing species. A summary is given in Section 5. 


2. Observations and Data Reduction 
2.1. General Observational Information 


The IRAM? Plateau de Bure Interferometer (PdBI) and 
National Radio Astronomy Observatory (NRAO) Very Large 
Array (VLA) observations in the 1.3 mm, 3.5 mm, and 1.3cm 
continuum toward the clumps G18.17, G18.21, G23.97N, 
G23.98, G23.44, G23.97S, G25.38, and G25.71 have been 
described and presented in Zhang et al. (2019). The NH2D 1,,- 
10; (at 85.926 GHz), NH; (1, 1), and (2, 2) observations have 
been presented in Zhang et al. (2020). Here the spectral 
observations and data, including HPCN ; HCPN, SO, and 
c!80, combined with the above data, are further described and 
exhibited. Additionally, in this paper, only the 3.5mm 
continuum data are used to trace the dense cores in the clumps. 
The 1.3cm continuum data, which mainly trace the free—free 
emission, are not used in this paper; 1.3 mm continuum data, 
which have relatively much smaller effective regions and lower 
sensitivity than 3.5mm continuum data, are also not used in 
this work. 


2.2. PdBI Observations 


The spectra in the PdBI observations were observed 
simultaneously with the continuum, but with separate correlator 
windows. In receiver 1, NH2D (85.926GHz) and HÜCN 
(86.340 GHz) were each tuned to a 40 MHz bandwidth with 
460 channels having a velocity resolution of 0.27 kms ! per 
channel; HCPN (86.055 GHz) and SO (86.094 GHz) were 
tuned to the same 80 MHz bandwidth with 230 channels 
having a velocity resolution of 1.09 kms! per channel. In 
receiver 2, C'?0 (219.560 GHz) was tuned to an 80 MHz 
bandwidth with 230 channels, which has a velocity resolution 
of 0.43 km s^! per channel. The rms noise of the spectra was 
about 23 mJy beam! for NH;D on the CD track and about 
12 mJy beam! for NH;D on the BCD track. In this work, we 
only used the observational data on the CD track. 

The IRAM software package GILDAS was utilized for 
spectral reduction. The region of a double primary beam was 
searched for cleaning components. No polygon was introduced 
to avoid any biased cleaning. The primary beam is 5875 and 
23"0 at 86.086 and 219.560 GHz, respectively. The data have 
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been corrected for primary beam attenuation. More details are 
presented in Zhang et al. (2019, 2020). 


3. Observational Results and Analysis 


For eight clumps G18.17, G18.21, G23.97N, G23.98, 
G23.44, G23.97S, G25.38, and G25.71, the maps for NH3, 
NHD, H"CN, HCPN, SO, C!%0, and the 3.5 mm continuum 
emission are shown in Appendix, while the maps for the kinetic 
temperature and column density can be found in Zhang et al. 
(2019, 2020). The clumps G18.17, G18.21, G23.97N, and 
G23.98 are in an early star formation stage or prestellar stage, 
while the clumps G23.44, G23.97S, G25.38, and G25.71 are in 
an evolved star formation stage or protostar stage (Zhang 
et al. 2019). We use the function scipy.stats.binned - 
statistic() in Python to count the pixel-by-pixel 
relationships between the tracers NH3, NH2D, H!’ CN, 
HCPN, SO, and C!8O and the 3.5 mm continuum emission, 
kinetic temperature, and column density. Each parameter space 
(from minimum to maximum) is divided into 20 bins. The 
minimum and the maximum of each parameter space can be 
seen in Figures 2-8, where we plotted some relationships of the 
tracers. In these figures above, the integrated intensities for the 
NH;, NH;D, H'°CN, HCPN, SO, and C'?O data are the 
averaged value in the given bin. Moreover, the 3.5mm 
continuum emission, kinetic temperature, and column density 
are relative values that have been normalized. The correlation 
for each tracer is shown below. 


3.1. NH; (1, 1) and (2, 2) 


We further present the NH; (1, 1) and (2, 2) spectral example 
shown in Figure 1 and integrated intensity maps displayed in 
Appendix. We mainly use the NH, data to show the 
distributions of the column density and kinetic temperature of 
dense cores. The calculations have been introduced in detail in 
Zhang et al. (2019, 2020). Furthermore, Zhang et al. (2020) 
presented the optical depth relations between NH;D and NH; 
(1, 1). We can see that most of the NH; (1, 1) cores have 
optical depth ranging from 1.5 to 9, indicating that these 
molecular clouds are often optically thick for NH; (1, 1). The 
optical depth of NH5D ranges from 0.1 to 9.7, most of which 
are more than 1, thus NH5D is usually optically thick in dense 
and cold cores. In Figures 2 and 3 we can see that the 
relationships between the NH3(1, 1), (2, 2) and the 3.5 mm 
continuum emission and column density are positively 
correlated for most of the sources. However, there is no 
obvious relationship between NH3(1, 1), (2, 2) and the kinetic 
temperature. In Figures Al and A2, we can see also that the 
NH; peak position is well consistent with that of the 3.5 mm 
emission, but the NH3 is much more extended than the 3.5 mm 
emission structure. However, the NH3(1, 1), (2, 2) distributions 
are more extended than the 3.5 mm emission distributions. 
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Figure 1. Example for spectra of H'^CN, HCPN, C!80, and SO extracted from the most massive 3.5 mm continuum cores (see the 3.5 mm continuum cores in Zhang 
et al. 2019). The green curves simply represent the line emission fit for each spectrum. 


[TTTTTT T[TTTTTTTTIETTTTTTTITTITTETTITTTTTITTTTTITE [T TTTTTTTTT[TTTITTTTT[TITTTTTTTITTTTTTTTITTT a| POT aa a E TTTTTTT[TTTTTTTTTTTTT] 
094r 3.5mm i 04r i 0.4 E Iun (1,07 NNB d 
E -=-~ G18.17 J E J E --*-- G18.17 E 
T E == GI821 na T me E =e G1821 4j @ E-*- GI821 k $ 
E 0.3 ---*- G23.97N Fae 4 B 03E -—- G2397N-] B 0.3 E---*- G23.97N FI Rf 
: p / wx a. : wee E E 
É -e-. G23.98 mo Ue m É æ- G23.98 É E -= 62398 oa 
Bl Se ee M rad S WON 000 824 $ o pna, E PRP 
2 —— 023975” e To s $ ZTO A ce eBo] $ be 9239 s a A Ay w* 
E 2p. 625.36 4 Mirat fetta E | m ope Gnaeet. ve a. 
= "m asi i / e-o E = "uu = E 4o tef ser 
iz 25 faethe’ A Ke = E Mis j 
T - i a rd we E eed $ J 
zl iy xti | E o arms = E ergas v. 
X 01 AA ee am ety "01 bE "X 01E EC ES 
- sar uo P N, x, E --e- m 
weg Sa N gt 3T T ee ean Sa E I ctu a eee. Wes cd 
ual unas i PU Ree aaa bara bia td 
. 0.6 0.8 $ 0.2 0.4 0.6 d 0.0 0.2 04 — 06 0.8 1.0 
Relative I5 sim Relative Tkin Relative Nun, 


Figure 2. Binned NH,(1, 1) integrated intensity distributions vs. related 3.5 mm continuum emission, kinetic temperature, and NH; column density distributions. 


3.2. NH5D at 85.926 GHz a few of the NH;D cores having kinetic temperatures above 

We further present the NHD integrated intensity maps in 20 K. It is likely that these cores with 220 K are just ahead of 
Appendix. Zhang et al. (2020) found that the NH2D cores have the deeply embedded cores with 720 K in line of sight. The 
a colder state than these continuum cores, and most of the derived kinetic temperatures may be contaminated by back- 
NH2D cores have a temperature between 13 and 17 K with only ground emission, and the NH3D is just tracing the envelope of 
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Figure 3. Binned NH3Q, 2) integrated intensity distributions vs. related 3.5 mm continuum emission, kinetic temperature, and NH3 column density distributions. 
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Figure 4. Binned NH;D integrated intensity distributions vs. related 3.5 mm continuum emission, kinetic temperature, and NH; column density distributions. 


the cores with 220 K. Therefore, the NHD lines may be 
depleted above 20 K. Comparing the column densities between 
NHD and continuum cores, we found that the NHD 
distributions are in a relatively diffuse condition compared to 
the continuum cores. Therefore, the NH5D line may be a good 
tracer of prestellar cores. In Figure 4 we can see that there is an 
anti-correlation between the NH5D and the kinetic temperature, 
but these are positively correlated for the sources G18.17, 
G23.44, and G18.21, while the relationship is not obvious for 
the other sources. In addition, the relationship between the 
NH3;D and the 3.5 mm continuum is also not obvious, but 
there is a weak positive correlation between the NH5D and the 
column density. The possible reasons for this are discussed in 
detail in Section 4.1. 


3.3. H?CN at 86.340 GHz 


H'3CN often traces dense and cold cores. Spectra of HCN 
have three hyperfine structure (HFS) lines (see Figure 1) 
(Csengeri et al. 2011; Padovani et al. 2011). The HPCN 
spectra of the most massive core among 3.5 mm continuum 
cores are shown in Figure 1. In Figure 5 we find that the 
relationships between the HCN integrated intensity and the 
related 3.5mm continuum emission, kinetic temperature, and 
NH; column density distributions are positively correlated for 


most sources. The integrated-intensity contours are super- 
imposed on the 3.5 mm emission image in Figure A5. We can 
see that the HCN emission is compact and surrounds the 
3.5 mm emission. It is likely that H'^CN often traces dense and 
cold cores for the 3.5 mm emission. 


3.4. HCN at 86.055 GHz 


The HCPN spectra of the most massive core among 3.5 mm 
continuum cores are shown in Figure 1. The integrated- 
intensity contours are superimposed on the 3.5 mm emission 
image in Figure A6. The HCPN distributions are well 
consistent with the 3.5 mm peaks. The cores within clumps 
G23.44, G23.97S, G25.38, and G25.71 are more evolved than 
those in clumps G18.17, G18.21, G23.97N, and G23.98, and 
the most massive cores have a relatively higher temperature 
condition than the surrounding cores (Zhang et al. 2019). 
Therefore, it is likely that the HCPN traces the evolved 
conditions, and the distributions of HCPN emissions are only 
concentrated close to warmer condition within each source. In 
Figure 6 we find that the relationships between the related 
3.5 mm continuum, NH; column density distributions, and the 
HCPN integrated intensity are positively correlated for most 
sources, but we cannot see a correlation between the kinetic 
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Figure 6. Binned HC'°N integrated intensity distributions vs. related 3.5 mm continuum emission, kinetic temperature, and NH3 column density distributions. 


temperature and the HCPN integrated intensity. This indicates 
that the HCPN is not a good temperature tracer. 


3.5. SO at 86.094 GHz 


The SO spectra of the most massive core among 3.5 mm 
continuum cores are shown only for sources G23.44, G23.97S, 
G25.38, and G25.71 in Figure 1, because we did not detect any 
SO emission at the other four sources. The integrated-intensity 
contours are superimposed on a 3.5 mm emission image in 
Figure A7. The SO distributions are well consistent with the 
3.5 mm peaks, except that the peak for G25.38 is offset from 
3.5mm emission. This molecule, SO, appears to be less 
confined to a small region than others. In Figure 7 we can see 
that there is an obvious positive correlation between the SO 
integrated intensity and the 3.5 mm continuum distributions. 
For the kinetic temperature and column density we do not see 
any obvious relationships. 


3.6. C50 at 219.560 GHz 


The C'8O spectra of the most massive core among 3.5 mm 
continuum cores are shown in Figure 1. c'8O emissions have 
been superimposed on the 3.5 mm continuum in Figure A8, 
where the primary beam of C50 is the dotted circle. Most of 


the C'*O emission is available at or near the 3.5 mm peaks 
because of the warmer conditions or more evolved condition 
there. It is quite rare that the NHD cores are located at the 
positions of the C'*0 emissions. Therefore, C'?O is susceptible 
to be severely depleted at low temperatures (<20 K) and high 
densities. In Figure 8 we can see a positive correlation between 
the C'^O integrated intensity and column density distributions 
only for source G23.44, while for other sources the relation- 
ships are not obvious. The main reason for this may be that the 
C'*0 is optically thick and has been frozen to the dust grain 
surface in the dense and cold condition (e.g., Zhang 
et al. 2017). Generally, we suggest that the C'^O could be 
used to study the diffuse (<10* cm ^) and warm (220 K) star 
formation regions. 


4. Discussion 
4.1. NHD Distribution 


NHD may form in different ways, for example, the 
deuterated ions trace the place of hydrogen in the ammonia 
in the gas phase, and in cold areas ($20 K), HD* reacts with 
NH; to generate NHD (e.g., Rodgers & Charnley 2001; 
Caselli et al. 2008; Caselli & Ceccarelli 2012; Ceccarelli 
et al. 2014). Li et al. (2024) presented NH2D mapping in a 
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Figure 8. Binned C'?^O integrated intensity distributions vs. related 3.5 mm continuum emission, kinetic temperature, and NH; column density distributions. 


sample of 24 late-stage massive star-forming regions using the 
IRAM 30 m, and found that NH5D emissions have a complex 
distribution. They suggest that the dense gas (or gas density) 
may not be an important physical parameter to affect the NHD 
enhancement in these targets, compared with other physical 
parameters (such as temperature). However, this may be biased 
due to their limited spatial resolution and sensitivity. 

In this work, we found that NH2D cores have quite narrow 
line widths, half of which are below 1.0kms ! (see Zhang 
et al. 2020), indicating a very small velocity dispersion for 
NH3D cores. Little NHD emission is detected toward the 
continuum cores, but we detected NH; lines for most of the 
cores (see Figure A4). Based on this, we can compare the 
velocity dispersion between NH5D and continuum cores (see 
Zhang et al. 2019, 2020). We found that the continuum cores 
have a larger velocity dispersion than the NHD cores. The 
thermal broadening of the spectra is very small at the rotational 
temperature below 20 K, so the turbulent conditions of the 
NH5D cores are relatively quiescent. 

For the eight samples in Figures Al and A3, we super- 
imposed NH5D and NH; contours on the 3.5 mm continuum. 
We found that the NH5D peaks are not associated with either a 
dust continuum or an NH; peak. For clumps G18.17, G18.21, 
G23.97N, and G23.98 in a prestellar stage, there are very weak 


infrared and millimeter emissions that were blended with the 
NH»5D distributions. For clumps G23.44, G23.97S, and G25.38 
in a protostar stage, the NHD distributions are extended and 
surround the 3.5 mm peaks. Only for clump G25.71 (in a 
protostar stage) is the NHD core located at the peak of the 
3.5 mm continuum. It is likely that the 3.5 mm continuum core 
is just located in the background of NH5D core along the line of 
sight. 

Based on the temperature distributions of the dense cores, we 
found that the suitable habitat (or living environment) for 
NH5D is 13-20 K, and the deuterium fractionations will reach a 
maximum at 16.5 K. For temperatures higher than 20 K, the 
activity of NH2D is likely to be inhibited. At temperatures 
below 13 K, NHD may also tend to be frozen out on dust 
grains. However, we should also consider that the rotational 
temperature tracer, NH3, may have been completely frozen 
onto the dust grain before NH5D. 


4.2. N-bearing Species 


Some N-bearing species are widely used as good tracers of 
dense gas and cold conditions for initial star formation, since 
they do not freeze out on dust grains in that condition (e.g., 
Suzuki et al. 2018; Awad & Shalabiea 2020; Peng et al. 2022). 
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However, different N-bearing species are slightly different in 
probing different physical conditions, relatively. In this work, 
we collect five N-bearing species: NH; (1, 1) and (2, 2), NH5D, 
HCN, and HCPN. From the integrated emissions in Figures 
A1—A3, A5, and A6, NH; (1, 1) and (2, 2) are more widely 
distributed than others, while HCPN lines are much weaker 
and appear only at the peak positions of these samples. 
Furthermore, we detect almost no HCPN emission in the 
direction of G18.17, G18.21, G23.97N, and G23.98, which are 
at the relatively early star formation stage (or prestellar stages), 
but we detect HCPN emission in G23.44, G23.97S, G25.38, 
and G25.71, which are in a protostar stage. Therefore, HCPN 
may be tracing these evolved protostars. NHD emissions are 
extended, but with an offset from the continuum and NH3 
peaks, so it can trace the earliest star-forming conditions among 
the five N-bearing species. H'^CN emissions are distributed in 
an extended way and agree well with the peak position of the 
3.5 mm continuum in all samples except G23.98. Therefore, 
HCN may cover a relatively extensive evolutionary process 
from early to evolved stage. 

Furthermore, Li et al. (2024) also found that NH2D and 
HÜCN have different distributions with asymmetrically and 
resolvably distributed spatial structures in 11 of 18 sources, but 
the others have no significant differences. They think that this 
is mainly because of their limited spatial resolution and 
sensitivity. Actually, in our high-spatial-resolution and high- 
sensitivity observations, the emission distributions between the 
NH;D and H'°CN are strongly different. This should result 
from their different habitats. The NH5D prefers to be located in 
colder conditions than the HPCN. 


5. Summary 


This work is a follow-up to Zhang et al. (2019, 2020), where 
we studied the gas dynamics, NH5D chemistry, and fragmenta- 
tions in eight massive precluster and protocluster clumps, 
including G18.17, G18.21, G23.97N, G23.98, G23.44, 
G23.97S, G25.38, and G25.71. The observational data are 
from the high spatial resolution interferometric observations 
with the PdBI and the VLA. In this work, we present PdBI and 
the VLA spectral data including NH3, NH»D, H?CN, HCPN, 
SO, and CO. We also study the possible habitats for these 
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tracers in the early stages of high-mass star formation. Some 
findings are following: 

The NH; cores are in good agreement with the 3.5 mm 
emission peak, but the NH4 is much more extended than the 
3.5 mm emission structure. The NH5D cores are mainly located 
in the temperature range of 13.0-20.0 K, and the NH3D lines 
may be strongly depleted above 20 K. Comparing the column 
densities between NHD and continuum cores, we found that 
the NH;D distributions are in a relatively diffuse condition. In 
these relatively early star formation conditions, we detect 
almost no HCPN emission in the direction of G18.17, G18.21, 
G23.97N, and G23.98, but we detect HC'°N emission in 
G23.44, G23.97S, G25.38, and G25.71, which are in a 
protostar stage. It is likely that HCP^N traces these evolved 
protostars. H'^CN emissions are distributed in an extended way 
and agree well with the peak position of the 3.5 mm continuum 
in all samples except G23.98. The SO distributions agree well 
with the 3.5 mm peaks for the evolved star formation stage, 
except that the peak for G25.38 is offset from the 3.5 mm 
emission. However, we did not detect any SO emission in the 
four sources of the prestellar stage. We also found that at low 
temperatures (<20K) and high densities (2109 em ?), the 
CO molecule rapidly freezes onto the dust grains and 
becomes severely depleted. Based on the above, CO is a 
poor tracer in conditions of very early star formation or 
prestellar cores, but CO could be used to study diffuse 
(<10* cm ^) and warm (220 K) star formation regions. 
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Figure A1. NH; (1, 1) integrated-intensity contours overlaid on 3.5 mm continuum emission with velocity range covering only the main line. The contour levels start 
at —3c in steps of 3c for NH; (1, 1) with c4. (s) = 12.5, 20.2, 9.1, 9.8, 15.8, 12.6, 12.9, 10.2 mJy beam"! km s~!. The synthesized beam size of each subfigure is 


indicated at the bottom-left corner. 
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Figure A2. NH; (2, 2) integrated-intensity contours overlaid on 3.5 mm continuum emission with velocity range covering only the main line. The contour levels start 
—3oc in steps of 3c for NH3 (2, 2) with o(4 (s) = 10.8, 13.5, 7.9, 7.2, 12.3, 8.2, 9.7, 8.3 mJy beam~'kms~!. The synthesized beam size of each subfigure is 


indicated at the bottom-left corner. 
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Figure A3. NH;D integrated-intensity contours overlaid on 3.5 mm continuum emission with velocity range covering all the six HFS lines. The contour levels start at 
—3oc in steps of 30 for NH5D with o@_) = 60.4, 60.9, 51.1, 53.2, 84.1, 48.3, 49.0, 50.4 mJy beam ! km s. The green numbers indicate the positions of extracted 
NH2D cores and condensations. The synthesized beam size of each subfigure is indicated at the bottom-left corner. 
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Figure A4. NH,D integrated-intensity contours overlaid on an NH; integrated-intensity image with velocity range covering all the six HFS lines. The contour levels 
start at —3c in steps of 3o for NH2D with oa) = 60.4, 60.9, 51.1, 53.2, 84.1, 48.3, 49.0, 50.4 mJy beam ! km s^ !. The synthesized beam size of each subfigure is 


indicated at the bottom-left corner. 
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Figure A5. H'°CN integrated-intensity contours overlaid on 3.5 mm continuum emission with velocity range covering all the three HFS lines. The contour levels start 
at —30 in steps of 3c for HCN with Ta-aa) = 63.3, 58.4, 37.6, 60.9, 97.6, 73.0, 84.8, 90.1 mJy beam"! km s~!. The synthesized beam size of each subfigure is 


indicated at the bottom-left corner. 
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Figure A6. HCPN integrated-intensity contours overlaid on 3.5 mm continuum emission with velocity range covering only the emission line. The contour levels start 
at —3oc in steps of 3c for HCN with Taa) = 88.5, 69.3, 73.8, 81.8, 61.4, 58.0, 55.1, 48.3 mJy beam"! km s~!. The synthesized beam size of each subfigure is 


indicated at the bottom-left corner. 
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Figure A7. SO integrated-intensity contours overlaid on 3.5 mm continuum emission with velocity range covering only the emission line. The contour levels start at 
—3o in steps of 30 for SO with e s; = 71.4, 61.4, 52.1, 46.7 mJy beam ! km s^. The synthesized beam size of each subfigure is indicated at the bottom-left 


corner. 
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Figure A8. C'^O integrated-intensity contours overlaid on 3.5 mm continuum (BCD track) with velocity range covering only the emission line. The contour levels 
start at —3c in steps of 30 for C!50 with Oaa) = 330.3, 151.1, 185.1, 204.7, 520.1, 214.6, 184.9, 198.8 mJy beam! km s~!. The synthesized beam size of each 


subfigure is indicated at the bottom-left corner. 
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